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We observe stationary random tessellations X — {'E.n}n>i in R'' through a convex sampling 
window W that expands unboundedly and we determine the total (k — l)-volume of those 
(fc — l)-dimensional manifold processes which are induced on the fc-facets of X (1 < fc < d — 1) 
by their intersections with the (d — l)-facets of independent and identically distributed motion- 
invariant tessellations X„ generated within each cell H„ of X. The cases of X being either a 
Poisson hyperplane tessellation or a random tessellation with weak dependences are treated 
separately. In both cases, however, we obtain that all of the total volumes measured in W are 
approximately normally distributed when W is sufficiently large. Structural formulae for mean 
values and asymptotic variances are derived and explicit numerical values are given for planar 
Poisson- Voronoi tessellations (PVTs) and Poisson line tessellations (PLTs). 
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tessellation; Poisson hyperplane process; Poisson- Voronoi tessellation; weakly dependent 
tessellation 

1. Introduction 

In this paper, we consider stationary random tessellations X = {S„}„>i of the d- 
dimensional Euclidean space IR'' with convex cells S„. We assume that within each cell 
S„ of the initial tessellation X, a further random tessellation X„ = {E!„£}£>i of is 
nested, that is, S„ is subdivided into cells S„ H 'B.nt,(. > 1, where the sequence of compo- 
nent tessellations (X„)„>i consists of independent copies of a generic motion-invariant 
tessellation Xq drawn independently of X . The assumption of motion-invariancc of Xq 
will play a crucial role in deriving explicit moment formulae. This type of iterated ran- 
dom tessellation is said to be an X/XQ-nesting in R'^. Having available only a single 
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observation of such an X/Xg-nesting in a presumably large, convex sampling window 
W , we are interested in the asymptotic behaviour of the random sums 

zf{W) = Y,^^^\W), l<k<d-l, (1.1) 

n>l 

where 'asymptotic' means that "f R'' and where the random measures 

^^i'H-)-'^fc-i(x('^"''nsWn-), i<fc<d-i, (1.2) 

act on the Borel sets of . The functional d';^\w) measures the (fc — l)-volume of 
the random subsets induced in W by the intersection of the motion-invariant manifold 
process X^~^^ of (d— l)-facets of X„ with the union of all fc-faces belonging to the 
boundary 9S„ of the nth cell of X (cf. Section 2.1 for details and precise definitions). 
Note that, by definition, the random measures i?J''^(.),4''^ (•),... in (1.2) are conditionally 
independent given the tessellation X. 

Our results supplement earlier central limits theorems (CLTs) for cumulative mea- 
sures of stationary ergodic tessellations modelling the total effect of random internal cell 
structures ([14]). Whereas, in the latter reference, the random measures corresponding to 
those in the sum (1.1) act on the interiors of the cells S„, the measures defined m 

(1.2) are concentrated on the cell boundaries 9E!„ of X. Hence, certain new effects arise 
due to the interactions between the stationary random manifold process IJri>i ^^^^ 
boundaries of X and the component tessellations {Xn)n>i- It turns out that there are 
considerable differences between X being a stationary Poisson hyperplane tessellation 
(PHT) and X satisfying certain weak dependence assumptions. In the first case, due to 
the overnormalization in the CLTs for Poisson hyperplane processes caused by inherent 
long-range dependences, (cf. [15]), the influence of Xq on the Gaussian limit distribu- 
tion is relatively weak. The other case seems to be somewhat more delicate because the 
asymptotic variance of the existing Gaussian limits are influenced by first and second 
order characteristics of both X and Xq. 

We present our derived CLTs in the general case of M"^, since this allows for a clearer 
and more transparent exposition. Clearly, however, the CLTs find their applications in 
the modelling of planar, but also spatial, networks as they occur, for example, in cell 
biology and telecommunications. Indeed, concentrating on the latter example, the prob- 
lem often arises of handling and modelling data that represent the geometrical struc- 
ture of the infrastructure system (e.g., main roads and side streets) that supports the 
technical telecommunications equipment. In recent years, stochastic-geometric modelling 
approaches have proven useful and are established domains of research today. In partic- 
ular, the Stochastic Subscriber Line Model (SSLM) has been developed as an integrated 
and easily extendable model for telecommunication access networks (cf. [25] and the 
references therein). 

The SSLM employs (iterated) random tessellations to describe the geometric network 
support. Having identified the best fitting model from a class of potentially suitable 
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(a) X = FLT (b) X = PVT 



Figure 1. Realizations of planar tessellations X, where the nested tessellation Xo is a PLT. 
The T-crossings are displayed as thick dots •. 

tessellations (cf. [8]), cost functionals and their distributions can be studied along the 
network geometry (cf. [9]). 

Assume that we use a planar X/Xg-nesting to model the geometric support. In the 
framework of our study, we observe a single value Z\ {W), which counts the number of 
T-crossings in a sampling region C induced by the intersection of the edges of the 
tessellation X„ with the edges of the nth cell S„ for n>l. Figure 1 shows two examples 
of this situation in differently shaped sampling windows W. In particular, Figure 1(a) 
shows a PLT/PLT- nesting through a ball of radius r > and centered at the origin, 
whereas, in Figure 1(b), we consider a PVT/PLT- nesting within a rectangular sampling 
window. 

The analysis of the aforementioned T-crossings, that is, of the connections between 
main roads and side streets, plays an important role in telecommunication modelling since 
these crossings are the entry points to the blockwise civil engineering of the local network. 
Let the type of the initial tessellation X and the type of the nested tessellation be 
known. Within a suitably large region W ^ the distribution of the number of T-crossings is 
then known through our results. Thus, the engineer is provided with useful information 
about the local network. For example, it is possible to deduce the dimensioning and 
capacity potential for each entry point in order to provide blockwise optimal connection 
quality to the subscribers, where one block comprises all of those subscribers who are 
situated in the cells formed by the main roads. 

In contrast to that, assume that we have, again in a suitably large region W , knowledge 
about local information, such as the type of Xq, and especially about the value Z{ {W) 
for the T-crossings. The expression on the left-hand side of (4.1), as well as the expression 
on the left-hand side of (5.7), can then be calculated and used to test for normality. 



CLTs for functionals on facets of random tessellations 



871 



Depending on the (unknown) type of X (representing the main road system) , we expect 
to reject the null hypothesis of normality either for the formula in (4.1) or for the formula 
in (5.7). This can provide, in the framework of model selection, a hint as to the structure 
of X before passing to more refined fitting procedures. 

The paper is organized as follows. In Section 2, we introduce basic notation and recall 
some relevant facts from stochastic geometry. Section 3 presents mean value relations 
and formulae for (asymptotic) variances. In Sections 4 and 5, we formulate and prove the 
announced CLTs for the different cases of initial tessellations X. Finally, in Section 6, we 
study some examples of weakly dependent tessellations and discuss possible extensions 
of our results. 



2. Preliminaries 

In this section, we introduce the basic notation and present a brief account of some rele- 
vant material on random tessellations and stochastic geometry in general. For a detailed 
and rigorous discussion of these topics, wc refer to the existing mathematical literature, 
in particular to [21, 22, 24, 27] and [30], which contain many further references, as well 
as numerous tessellation models with applications to various fields. 

Throughout, let (f2, (T(f2), P) be a common probability space on which all random ob- 
jects occurring in the present paper will be defined. Further, let {x,y) = J^t^i^kVk de- 
note the scalar product of the coordinate vectors x = (xi, . . . , x^)^ and y = {yi, . . . , yd)^ 
in . By means of the Euclidean norm || • || = \J {■, •), wc define the closed ball 
B'^ — {x ^'R'^: \\x\\ < r} with radius r > centered at the origin and the unit sphere 
S''""^ = {x G R'':||a;|| = 1} in M'', respectively. Remember that each afhne {d — 1)- 
dimensional subspace H of M.'^, called hyperplane in Mf^ in the sequel, admits a parameter 
representation H{p,v) = {x GMf^ : {x,v) =p}. Here, p S denotes the signed perpen- 
dicular distance of H from the origin and v G S"^^ = {{xi, ■ ■ ■ , Xd)'^ S S''"^ : > 0} is 
the directional vector belonging to the upper unit hemisphere. Further, let J^fc(-) de- 
note the Lebesguc or fc- volume measure in MJ^ for k = 0, . . . ,d, where we can also just 
write i^d{') = I • I- The /c-dimensional Lebesgue measure will also be used instead of 
the /c-dimensional Hausdorff measure on (afhne) fc-dimensional subspaces in M.'^ for any 
fc = 0, . . . ,d— 1. As usual, vq^-) coincides with the counting measure, that is, vq^B) = jj=B. 
For brevity, put 



The families of all non-empty closed, compact and compact convex sets in R'' arc denoted 
by T'^, JC'^ and C^, respectively. Note that M{S) stands for the a-algebra of Borel sets in 
the metric space S. 




.d/2 
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2.1. Random tessellations and random nestings 

In this section, wc sketch out the mathematically rigorous approach to random tessel- 
lations, as used in stochastic geometry, and we recall some basic facts, where referring 
[21, 22, 24, 27] and [30] for a systematic study of these topics. 

A tessellation of is a countable family r = {C„}„>i of convex bodies C„ e C'^ 
such that intC„ ^ for all n, intC„ n intCm = for all m, [J^^i C„ = R'^ and 
Sn>i l{c„nA'#0} < for any K G /C^. Notice that the sets C„, called the cells of r, 
are necessarily polytopes in R''. The family of all tessellations in R'^ is denoted by T. A 
random tessellation X ~ {S„}„>i in M.'^ is a sequence of random convex bodies S„ such 
that P(X e T) = 1. 

Note that a (stationary) random tessellation X can also be modelled as a (stationary) 
marked point process J2n>i ^[a{E„),E°] > where a : ^ M'' is a i3(^^)-measurable mapping 
such that a(C) S C and a{C + x) = a{C) + x for any C G C'^ and x S M'^, and where 
S° = S„ — a(E!„) is the centered cell corresponding to S„ which contains the origin. The 
point a(C) is called the associated point of C and is usually chosen to be the centroid or 
lexicographically smallest point of C. 

Suppose that the stationary marked point process X]?i>i '^[q(s„),ho] has positive and 
finite intensity 7 ~ E#{ri : a(S„) G [0, 1)''}. By V^, we denote the set of all compact and 
convex d-polytopes whose associated point is located at the origin. The corresponding 
Palm mark distribution P'^ of X is then given by 



The notion of a typical cell of X refers to a random polytope S* : ^ whose distri- 
bution coincides with P*^. Since the cells S„ are space-filling and non-overlapping (up to 
a null set), we have the mean value relationship 



that is, the cell intensity 7 equals the reciprocal of E|S*|. 

A (deterministic) iterated tessellation r = {C„f n C„ : int Cni H int C„ 7^ 0} in con- 
sists of an initial tessellation t = {C„}n>i in M'' and a sequence (t„)„>i of component 
tessellations t„ = {Cni}e>i- In order to define a random iterated tessellation, we proceed 
along the lines of [19]. Let S be a random convex body in E'^ with P-a.s. non-empty in- 
terior and let X = {S„}„>i be a random tessellation in M.'^. Then, the counting measure 



Y{- I S) defined by Y{B \ S) = E„>i '5s„nH(-B)l|i„tH„nintB^0} for B e B{r^) is a point 



process in describing a random tessellation of S. 

Furthermore, if X = {S„}„>i is an arbitrary random tessellation in and if Xn = 
{S„^}£>i, n= 1,2..., are independent copies of a generic random tessellation Xq in 
M'^ drawn independent of X, then the random counting measure Y{B) = Yn{B \ S„), 
where Y„(B \ S„) = J2e>i fe„«nE:„(-B)l{i„tH„jnintH„5:^0} for B e B{T^), is called the point- 
process representation of an iterated random tessellation (briefiy , an X/ X^-nesting) in 



P\B) = 7'^E#{n : a(S„) e [0, 1)^ Si e B}, 



BeBiT'^)nvl 



(2.1) 




(2.2) 
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R with initial tessellation X and component tessellations Xi,X2, .... Clearly, the point 
process Y is stationary (and isotropic), provided that both the initial tessellation X and 
the generic component tessellation Xq are stationary (and isotropic). Moreover, Y is 
ergodic if X possesses this property. 

Each stationary (motion-invariant) random tessellation X = {S„}„>i in Mf^ induces d 
stationary (motion-invariant) random lower-dimensional manifold processes X'^''\ called 
k-facet processes of X for fc = 0, — 1. For example, X^^^ is the point process of 

vertices and X^^^ is the line segment process of edges of X . 

To be precise, X^'^^ is defined to be the union of all of the k-facets of X, whereas si*^'' 
denotes the union of all k- faces of its nth cell S„ . Here, the fc- facets of X are fc-polytopes 
in R'^ which arise from a finite intersection of neighbouring cells of X. The {d— l)-faces 
of S„ are (d— l)-polytopes in the boundary 9S„ and fc-faces arc defined recursively for 
fc = 0, . . . , d — 2 as fc-polytopcs in the relative boundaries of the (fc + l)-faccs. Note that 
the set of all fc-faces may differ from the set of fc-facets and that, for example in [27], 
Chapter 6, X^''^ is used slightly differently to denote the point process of fc-facets. 

A random tessellation X = {S„}„>i in M'' is said to be normal (or ordinary) if P- 
a.s. every fc-facet of X lies in the boundaries of exactly d — fc + 1 cells, k ~ 0, . . . ,d ~ 1. 
Many-real life tessellations in and possess this property, which motivates the 
term 'normal'. There are important classes of stationary tessellations in R'' whose cells 
are constructed (realizationwise) , according to specific geometric rules, from the atoms 
of a stationary point process in R''. Among them are Voronoi and Laguerre tessella- 
tions (see, e.g., [24] for details), which turn out to be normal if the generating point 
process is Poisson; sec [21]. It seems that this fact continues to hold for a large class 
of (even instationary) generating point processes which are mixing in a certain sense 
and/or whose higher-order moment measures possess Lebesgue densities. In [13], it is 
shown that Voronoi tessellations in R'^ are normal if the {d + 2)th-order product density 
of the generating stationary point process exists. Computations of second-order charac- 
teristics of spatial Poisson- Voronoi tessellations can be found in [11]. Finally, it should 
be mentioned that there are more general definitions of tessellations (cf., e.g., [32]), al- 
lowing for the rigorous treatment of random tessellations which do not necessarily consist 
of only convex cells. Without doubt, the most prominent example is the Johnson-Mehl 
tessellation; see also [24] for details. For this model, CLTs have been proved, based on 
a-mixing conditions derived from the generating (Poisson) point process (cf. [4] and [5]). 

2.2. Stationary Poisson hyperplane tessellations 

Let 4' = X)i>i ^[Pi-Vi] be a stationary and independently marked Poisson point process on 
the real line R^ with intensity A and mark distribution 8 on the mark space S!^~^; see [6]. 
By means of the parameter representation H{p, v), {p, v) <eS} x S*^"^, of a hyperplane in 
R'^, we may represent a Poisson hyperplane process (PHT) <!> (defined in [27] as a point 
process on the space of afhne {d— l)-dimensional subspaces in R'') with intensity A and 
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(spherical) orientation distribution & by 

I'-EW-.v.)- (2-3) 
j>i 

The Poisson hyperplane process $ given in (2.3) is said to be non- degenerate if 
&{H{0,v) n §+~^) < 1 for any v G S'fr^ In this case, (2.3) induces stationary k-fiat 
processes <f>fc for fc = 0, 1, . . . , c? — 1 whose countable support consists of the affine k- 
dimensional subspaces {k- intersection flats) H{Pi^ , K: ^ ) H • n H{Pi^_^ , Vi^ ,,) for pairwise 
distinct indices ii, . . . ,id-k > 1- The union of these fc-flats coincides with the fc- facet pro- 
cess X^''^ of the corresponding stationary PHT X = {S„}„>i generated by (2.3). The 
cells S„, n>l, are bounded d-polytopes (P-a.s.) if and only if $ is non-degenerate; see 
[27], Chapter 6. Furthermore, this property implies that the stationary fc-volume measure 
'&k,d{-) associated with $fe (resp. X'^'^^) and defined by 



^kAB) = ,^ ,,, y* ^ '^M n H{P^, ,V^^)r\B] for bounded B e B(R^), 

(d — fc)! ^^n,...,4d-fc>i \ _[ J 

_ ^ (2.4) 

where ^ denotes summation over pairwise distinct indices, has positive intensity 



Afc.d^E^fc,rf([0,l)'^)^ EgkAQo,Vo) for = 0, 1, . . . , d - 1. (2.5) 

Here, the function {p,v) g]^\p,v) is defined on [—1, 1] x "B*^^ by 

(d-k-l \ 
Q H{Q,,V,)nH{p,v)nBfV (2.6) 

where {Qi, Vi), i = 0, 1, . . . , d — 1, are i.i.d. random vectors with independent components. 
Note that the generic random variable Qo is uniformly distributed on [—1,1] and the 
generic random vector Vq has the orientation distribution Q; see [15]. 

It is well known from convex geometry (see [26]) Chapter 3.5, that the probability 
measure Q on S!|~^ determines a unique, centrally symmetric convex body Zq, called 
the associated zonoid, which is given by 

h{Ze,u)= [ |(w,w)|e(dt;) foruSR'', 

where h{K,u) = Taa.XxeK{u,x) denotes the support function of an arbitrary K e C^. 

In [16], the following closed-form expression of gk,d{p,v) in terms of Zq has been 
derived 

9kAP,v) = ^^',',:,'_^!"'^-^ l-P^)''^-^^/^lhMl(rt^i^;!\(^e), (2.7) 
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where denotes the image of K G C'^ under orthogonal projection onto H{0,v) and 
yjd stands for the intrinsic j-vohimc of K <E C^_i- Using the relationship 

jVj''\Ze)= [ vj^~'\z^Mdv) forj = l,...,d 

(cf. [26, 31], Chapter 3.5 and [16]) combined with (2.5) and (2.7) yields that XkM = A''"'^ 
X vj'^i^{ZQ) for /c = 0, 1, . . . , d — 1, which has already been stated in [20], Chapter 6. 

The stationary Poisson hyperplanc process <f> given in (2.3) is isotropic (and hence 
motion-invariant) if and only if is the uniform distribution, which means that Zq = 
^^^Bf. This, in turn, leads to the explicit formula 



kj Kk V dKd 



for fc = 0,l,...,(i- 1. (2.8) 



We are now in a position to formulate a CLT for the total fc- volume -dk.diBf) of the 
support of the fc-flat process $fc contained in the ball B'^. This result has been proven in 
[15], even in a multidimensional version. 

Theorem 2.1. Let $ = X]i>i ^6 o, stationary, non- degenerate Poisson hyper- 
plane process with orientation distribution O on Ei^^ and intensity A > 0. Then, 

^k,d{B'',)-\kABt\ d ,ru. 2 



^d|l-l/2d 



N{Q,ctI,) for fc = 0,l,...,d-l, (2.9) 



where 



\s.Mt\Bi)) (2A)2''-2*^-i 
aid = lim d\2-i/d = 2-i/d %M(Qo,^o), (2.10) 

with gk.diViV) and {Qq,Vq) defined by (2.6). If, additionally, $ is isotropic, that is, O is 
the uniform distribution on S'jr^, then Xk.d is given by (2.8) and cr^ ^ takes the explicit 
form 



Note that, even in the anisotropic case, we have 



"'-^■^ = '(2^-1)!.^/- ^'-''^ 

that is, (J^_i d coincides with the left-hand side of (2.11) for k = d — 1. This is due to 
the fact that the {d— l)-volume of the hyperplanes H{Pi,Vi) within the ball B'^ does 
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not depend on Vi and so the distribution of '0k,d{Bg) is independent of the orientation 
distribution &. Furthermore, we mention that in [16], Theorem 2.1 could be extended 
to non-spherical convex sampling windows Wg = gWi (cf. Section 5 below). However, in 
this case, the formulae (2.7) and (2.11) depend on Wi and are less explicit. 



3. First- and second-order moment formulae 

Let X = {S„}„>i be a stationary random tessellation of R'^ and let Xq be a motion- 
invariant tessellation independent of X. We consider an X/Xo-nesting in M'*, as in Sec- 
tion 2.1, observed within a convex sampling window W . In order to calculate expectation 
and variance of the random variables zjf"^ in (1.1) for fc = l,...,(i— 1, we need two in- 
tensity values. 

First, wc consider Xq^''^\ the intensity of the stationary (fc — l)-dimcnsional manifold 
process Xq'^ n i generated by the intersection of the {d — l)-facet process Xq'^ 
with an arbitrary fc-flat L in R''. Since Xq'^ is motion-invariant by the assumption of 
motion- invariancc of Xq, wc may identify L with M*"' so that Aq'^''*'' can be defined by 

A[,'''''^==Ez/fe_i(X^''"^^nR'^n[0,l)''% l<fc<d-l. (3.1) 

By using quite general stereological relationships derived in [18], we may express Aq'^''''' 
by the (full-dimensional) intensity Ag'^''^^ = Ei'ii-i{Xq'^~^^ n [0, 1)'^) of the manifold process 
of {d— l)-facets Xq'* through 

.{k.d)_(d).{d,d) (d)_ 7((fc+l)/2)7(rf/2) , , , , ,oo^ 

Ao -c, Ao withc, l<fc<d-l. (3.2) 

Further, let 

//i'')=Ez.fc(xWn[0,l)'*) (3.3) 

denote the intensity of the stationary fc-facet process X^*^^ associated with X. To avoid 
rather involved formulae, in particular, for the variance of zlf\w), we impose an addi- 
tional condition on the tessellation X = {5„}„>i. 

Condition F. For fc = 1, . . . , d~ 1, assume that there exists a non-random integer m^^'^ > 
1 such that 



-a.s. 



for any W e C'^ with \W\ > 0. 

Condition F means that, for k = l,...,d— 1, each fc-facet of X lies in a constant 
number m'jf'' of fc- faces of cells E!„ , rt > 1 . Obviously, Condition F is satisfied for any 
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(2) 

planar tessellation with m\ ' =2. For d > 3, by the very definition, any normal tessellation 

X obeys Condition F with mlf"^ = d — fc + 1 , and a non-degenerate stationary Poisson 

hyperplane tessellation X is easily seen to satisfy Condition F with m^jf'^ = 2''"'^, see [27]. 
Note, however, that Poisson-Delaunay tessellations do not satisfy Condition F for d>3. 

Lemma 3.1. Consider an X/X^-nesting in M'' with stationary initial tessellation 
X = {E!„}„>i and motion-invariant component tessellation Xq. Assume that X satis- 
fies Condition F and that < 7^"'^ =E|S*| < 00 (cf. (2.2)). If ^J■]^^ < 00 or, equivalently, 
Ei/fe(S*'-'^'') < 00, and Aq'^'''' < 00 for any k ~ 1, . . . ,d ~ 1, then 

Ezlf\w)^X^^^''^'>m[^^n[^'^\W\ foranyWeC'a, A: = 1, . . . , d - 1. (3.4) 
Moreover, if additionally, 

Ei/^(X('=)n[0,l)'')<cx) and [ E!.^_i(X^''"^' n C('^')P"(dC) < cx), (3.5) 

then for W €z C'^ and k ~ 1, . . . ,d— 1, 

Yarizlf\w)) = -/ [ [ Yar{iyk-iiXlj~^ (1 C^'''^ n {W - x))) dxP°{dC) 

(3.6) 

+ {\'i'''^mff\a.v{vk{X^^^r\W)). 

Proof. Let fc€ {!,... ,d— 1} be fixed and let Ex (•) denote the conditional expectation 
E(-|X) given the tessellation X — {S„}„>i. Hence, we may rewrite the expectation of 

zjf\w) introduced in (1.1) as 

n>l 

Owing to the motion-invariancc of Xq"^ ^\ wc get, together with (3.1), that 

Exiyk^x^,"-'^ n si^) nw) = xi^'^^^M^'n^ n w) 

for any cell S„. In view of Condition F, wc may proceed by writing that 

ExZ^^\w) = A^,"^'') ^^(^n ' ^ ^) = ^i'^'^rni^^'MX^'^ n W). (3.7) 

n>l 

Combined with the stationarity of X^''\ this gives EukiX'-''^ nW) = ^l[^^\W\, which, in 
turn, proves (3.4). Recall that by using the notion of the typical cell (cf. (2.2)), we have 

Ej2'^ki^il:^(^W)=-fE,.,{E*<^'^^)\W\{=m^^hkjW\), 

n>l 
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which establishes the relationship Evki'^*''^''^) = m[^^ ^kM^^*]- To verify (3.6), we start 
with the well-known identity 

YaT{zlf\w))=E{VaTxZlf\w))+yaTiExZi''\w)), (3.8) 

where Varx(-) denotes the conditional variance Var(-|X) given X. Since, conditional 
on the tessellation X = {S„}„>i the random measures i?^'"-' (•), tJj'^'' (•),.. . in (1.2) are 
stochastically independent, we obtain 

Varx zi^'Hw) = ^ Ya.vx{:^k-i{X^o'''^ n S^f) n W)). (3.9) 

n>l 

With S„ = S° + q;(S„), we may apply the refined Campbell theorem to the stationary 
marked point process J2n>i '^[q(h„).e:o] (cf. [6] or [21]), where we find, together with (2.1), 
that 



E{YaTxZlf\w))=-/ [ [ YaT{vk-i{X^o'^~^'^ n{C^''^ +x)nW))P"{dC)dx 

jR<i J-p° 

(3.10) 

= 7/ / Ymiuk-iiX^o'^^^'^ nC^''^ n{W ~x)))dxP°{dC). 

J-p° JR'^ 

Here, we used (3.5), the stationarity of Xq"^ ^\ Fubini's theorem and the fact that 

j^k-iiiB + x) r^w)^ vk-i{B r^{w-x)) 

for any bounded B e S(R'^^^) and x E M''. The existence of the inner Lebesgue integral 
in the second line of (3.10) is also seen by applying Fubini's theorem and the second 
condition of (3.5), that is, 

/ EvI_^{x^q'^-^^ nc'^''^ n{w -x))dx 

= e/ / \iW-u)niW~v)\iyk-iidu)iyk-i{dv) 

<Ei/2_^(X^'^-i)nC('^))|Ty|. 

From (3.7), combined with the first condition in (3.5), it is immediately seen that the 
second term on the right-hand side of (3.8) is finite and takes the form 

Yai{Ex Zlf\w)) = {Xll'^''^mi^^ fYi^r{iyk{X^'''> n W)). 
The latter equality, together with (3.10) and (3.8), confirms the validity of (3.6). □ 
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The second condition of (3.5) imposes restrictions on both the initial and the compo- 
nent tessellation. Note that this condition is fulfilled if 

Ei^^_i(X^'*"^^nR''n [0,1)'') <oo and EN'^ {E* ) (1 + D{E*)f'" < oo (3.11) 

for fc = 1, . . . , d— 1, where Nk{C) and D[C) = sup{||a: — : a;,?/ e C} denote the number of 
fc-faces and the diameter of the d-polytope C GV^, respectively. To see that (3.11) implies 
the second condition in (3.5), we write C'^'^^ as union of the fc-faces \ / = 1 , . . . , iVfc (C) , 
and use the motion invariancc of Xq'' to obtain the estimate 

E^ti(x(''-i)ncW)<E( J2 ^fe-i(^r''nc('M 

< iVfe2(C)Ei/^_i(x^''"^^ n M'^ n [0, D{cyf) 

< Nl{C){i + D{C)f^Eul_^{x^'^~^^ nrn'^n [o, 1)'=) 

for k ~ 1, . . . ,d — 1. 

Note that if i?l'°^(-) acts on the interior of E!„, as supposed in [14], then the conditional 
expectation Ex'&rt^W) is a constant multiple of |S„ n W\ and therefore 'ExZjf\w) = 
J2n>i'^x''^7i\w) is proportional to In this case, the second variance term on the 

right-hand side of (3.8) vanishes. Due to this fact, the formula for the variance 

given in [14] is relatively simple and the ergodicity of the initial tessellation X suffices to 

prove asymptotic normality of zjf\w). 

4. CLTs for manifold processes on facets of a 
stationary PHT 

In this section, we consider the random measures given in (1.2) for n > 1 whose 

support lies in the cell boundaries more precisely in the fc- facets (for fc'=:l,...,d— 1) 
of a non-degenerate stationary PHT X = {S„}„>i in K''. For the sake of simplicity, we 
assume in this section that the sampling window W is the d-dimensional ball centered 
at the origin and with radius g> 0. For more general expanding sampling windows, we 
refer to the comment at the end of Section 2.2. By Theorem 2.1, combined with the 
geometric and probabilistic properties of PHTs, we prove a CLT (as g oo) for the total 
(fc — l)-volume of the sets (contained in B^) arising from the intersection of the fc-faces 
of S„ with the {d— l)-facets of the component tessellations X„ for n > 1. 

Theorem 4.1. Let X = {S„}„>i be the stationary PHT in R'^ generated by a station- 
ary non- degenerate Poisson hyperplane process $ given in (2.3) with orientation dis- 
tribution and intensity A > 0. Furthermore, let Xq be a motion-invariant random 
tessellation in having the intensities X^q''''^ > (cf. (3.1) and (3.2), resp.) Assume 
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that Ev'^_^{Xq'^^^'^ n R*"' n (0, 1]*^) < oo for fc = 1, . . . , d — 1 and that the corresponding 
X/Xo-nesting is observed through the spherical sampling window B^. Then, 

|ij^^|i-i/(2^) „-^-^(0'^m) fork=l,...,d^l, (4.1) 

where 

Aifc = 2 Aq Afc,d, (Tfc^d = Jnn = (2 Aq ) a^^^ (4-2) 

anc? Afc^rf, (T^ and X^''''^ are defined by (2.5), (2.10) and (3.1), respectively. 

Proof. Wc first recall that in case of a stationary PHT X, we have that m'"^^ = 2'*^'^ 
and that the intensity (3.3) of the fc-facet process X^'^^ coincides with the intensity (2.5) 
of the fc-flat process induced by (2.3), that is, we have that ^^'^^ = \k,d- Hence, the 
formulae for the intensities 'jif^ = 'EZ^\[Q,lY) of Z^\-) follow from (3.4), as stated in 
(4.2). Next, we rewrite the mean zero random variable Z^'''\B'^) — /I^'^^ji?^! as 

- r^f\Bt\ = 5f ) + 2'^-^-A('^^'^)r(^-), (4.3) 

where 

5f ^ = zi'\Bt) W\^\Bt) - z\^\Bi) 2'^-^\t''>.,{X^^^ n Bi) 

and 

) = ^,(xw n i?,^) - t,f\B''^\ = §kABt) - \k,d\Bt\. 
From Theorem 2.1, wc obtain that 



(fc) 



e 



N{Q,al,) forfc = l,...,d-l. (4.4) 



|^d|l-l/(2d) „_oc 

By means of Slutsky's theorem (cf., e.g., [17]), the proof of the CLT (4.1) is complete 
whenever 



Sq 



|^d|l-l/(2d) g^oo 

In view of Chebychev's inequality, we need only to prove that 



(fc) 

P — >Q. (4.5) 



=<^"'' ^0, (4,6) 
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which includes first of all to ensure that ¥.{sf^f < oo. Since E(S'f = E Varx(^f 
wc obtain, in analogy to the proof of Lemma 3.1 and by taking into account (3.11), that 

E(5f)2^7/ / VarK_i(Xo^-inC('=)n(i?;;-a;)))dxP"(dC) 



< 



7 / E^ti(x(^-^)ncW)p°(dC)|i?,^| 



■po 



< jEul_,ixt'^ nR"^ n [0, l)'=)EiV|(S*)(l + DiEnflB^.l 

Using the distributional properties of the typical cell S* of a stationary PHT. in par- 
ticular that has an exponential moment (cf. [2]), we find that 

EN^{E*){l + D{E*)f'' <oo for fc=l,...,d-l, 

which immediately confirms (4.6) for any d>2 (cf. also [7] and [28]). Finally, using the 
formula (3.6) for the variance of zlf\B'^), together with the limiting relations (4.6) and 
(2.10), we find that 

~2 _ 1- Var(4''^(gg)) _ fr,d-kAk,d).2 2 

This completes the proof of Theorem 4.1. □ 

5. CLTs for manifold processes induced on the facets 
of a stationary weakly dependent tessellation 

Throughout this section, we consider a stationary X/Xo-nesting which can be observed 
through an expanding family of convex sampling windows Wg with shape Wg = gWi 
for g > 0, where Wi S C'^ contains a ball and is itself contained in a ball, that is, 
B!^ C Wi C Bf^ for some < r < R < oo. We assume that the stationary initial tessellation 
X = {S„}„>i is ergodic (cf. [6, 24, 27]) and possesses, in contrast to Poisson hyperplane 
tessellations, further weak dependence properties. The latter properties ensure asymp- 
totic normality of the total A:- volume of the fc- facets in a large sampling window Wg. 
More precisely, we impose on X the following condition. 

Condition G. For k = 1, . . . ,d~ 1, assume that there exists a real number t| ^ > such 
that 

Var(z.fc(xWnW^e)) ^ 
\Wo\ e^^^^^ 



MX^'^^r)Wg)-f^f\Wg, „ 

„-^-^(0. 
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In analogy to Section 4, we shall prove that the centered and normalized cumulative 
functional (1.1) on Wg, that is, 

^^"'T^.^/^"^'' with,f =EZr([0,l)^), (5.1) 

converges in distribution to a Gaussian random variable Af{0, t| j) for k = — 1, 

as Wg t K''. Here, ^ denotes the asymptotic variance of the random variable (5.1) as 
g —f oo, that is 



^^^^^ Var(zf (Ty,)) 



^W-;™ ,.V^j " for/c = l,...,d-l, (5.2) 



the existence of which is shown in the subsequent lemma. 



Lemma 5.1. Let there be given an X / X^-nesting in with stationary (not necessarily 
ergodic) initial tessellation X and motion-invariant component tessellation Xq satisfying 
the assumptions of Lemma 3. 1 such that the asymptotic variance t| ^ in the first part 
of Condition G exists. Then, the asymptotic variance ^ in {5.2) exists and takes the 
form 

rla = {ri'^'^? + (A^mi'^') /or A = 1, . . . , d - 1, (5.3) 

where 

(rf'')2-7/ Var(^,_i(x(''-i)nCW))P°(dC). (5.4) 



Proof. The proof of (5.3) is based on the representation of the variance oiZl^^Wg) for 
fc = 1, . . . , c? — 1 given in Lemma 3.1. From (3.6) and the first part of Condition G, it is 
easily seen that (5.3) holds if and only if the limit 

EjY^r x{zi'\Wg)) 
= lim / / Var(i.fc_i(X^''-')nC('=)n(W^,-x)))d.TP"(dC) 

g^co \Wg\ J-pO J^d 

exists and equals (tq*'''''')^, as defined in (5.4). To show this, we apply the same arguments 
as those already used in the proof of Lemma 3.1 to derive the estimate 

E{V&ixZi''\w))<-f\W\ [ Ez.2_^(^^d-i)^^(fc)^_po^^^^^ 
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By multiple application of Fubini's theorem, we arrive at 



E / -^^^ -— ^dx P°(dC) 



Note that, in view of limg^oo \ (Wg - u) n (Wg - v)\/\Wg\ = 1 and \(Wg-u)n{Wg-v)\ < 
\Wg\, together with (3.5), we may apply Lebesgue's dominated convergence theorem. 
Likewise, we obtain that 



(E^fc_i(x<'^-^>ncWn(VF,-.T)))^ 



da:P"(dC) 



Jv^JR" \^q\ 

V(M)^2 /• f f \{Wg-U)n{Wg-V)\^^ ,,^.,^0, 



= iK"y / / " ^ ' ^^v ,(dz.K(d«)p"(dc) 

^(A^'^'^Z t^,'(CW)P°(dC)= / (Ei.fc_i(X^''-'^nCW))'p"(dC), 

which completes the proof of Lemma 5.1. □ 

Under Condition F, we may decompose the normalized cumulative functionals given 
in (5.1), in analogy to (4.3), as 

^'"";^f"^' -^f'^<Vrvy', (5.5) 



where 



r(k)_ Zi'\Wg)-E^{Zlf\Wg)) ^^,,,_ MX^'^nWg)-X[,''-'^\Wg\ _ 



\Wg\^l^ '^"^ " |M^,|l/2 



Notice the fact that, for fc = 1, . . . , c? — 1 and any fixed > 0, the random variables 
Uf^ and yj''^ are uncorrelated. Obviously, Vg^^ is a (measurable) function of X and 
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IEx(C/f = (P-a.s), so 

E([/f Vf')) = E(Ex([/f ))) = E(Ex(C/f = 0. 

The following theorem states that, for any = 1, the two-dimensional vector 

{Uq^\ Vq^^)^ converges in distribution to a mean zero Gaussian vector with independent 
components as oo. This, in turn, implies the desired asymptotic normality of (5.1). 

Theorem 5.1. Consider an X/Xo-nesting in Mf^ observed through the increasing family 
of windows Wg with motion invariant component tessellation Xq and stationary ergodic 
initial tessellation X = {S„}„>i satisfying E_D'^(S*) < oo, as well as (3.5), Conditions 
F and G Then, 



T, 





2 

kA 



for k = 1, . 



,d-l. 



(5.6) 



In particular, this implies that 



id) I 



m,rla) fork^l,. 



,d- 1, 



(5.7) 



where 



Ad) 



Ak,d)(d)^(d) 



'k.d - 



(rf )' + (A[,'''^mi'')^2^2 



k.d- 



Proof. We employ the method of characteristic functions (cf., e.g., [17] for details). 
Hence, we must show that the characteristic function fg{s,t) of the random vector 

(C/f\v^i''V defined by 

/,(s,i)-Eexp{isC/W+itFW} 

converges to the characteristic function of the Gaussian random vector that occurs as 
limit in (5.6), that is. 



IQK-^,") ' CXp 



-Ti. 



k.d 



for all s,teR\ 



For this, we introduce the decomposition fg{s,t) fg^\s,t), where 



fj,'\s,t)=E 



Ex( exp{ist/f' }-exp 



YaTxizi'\Wg))}]eM^tV^'>} 



2W„ 



Ak)^ 



exp 



2\Wo 



■YaTx{zi'\Wg)) 



exp 



4r(-r^)' exp{iiFf)} 
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/f (s,t) = exp 



In view of Condition G, the continuity theorem for (one-dimensional) characteristic 
functions yields that 

hm /(3)(s,t)^exp(-^(T('='''V| lim Eexp{itV^f)} 



exp 



for all e M^. Hence, it remains to prove f^\s,t) - 
subsequently show that 



'k,d 



^ Q — ^OO 



for I = 1,2. For this, we 



W„\ ^— ' g^oo 



(5.8) 



where the first equality follows from (3.9) and the almost sure convergence of the argued 
sum in (5.8) can be argued by some modified ergodic theorem for random tessellations 
(cf. Theorem 4.1 in [14]), which states that 



\W, 



(5.9) 



for any yS(J^^)-measurable, translation invariant set function g{-) defined on sets of the 
form C n W, where C is a d-polytope and W G C'^^, and satisfying the monotonicity 
property g{Cr\ W) < g{Cr\ W) for W C W' . It is easily verified, by checking the proof of 
Theorem 4.1 in [14], that these restrictions imposed on g(-), together with EZ3'^(S*) < oo 
and Eg{E*) < oo, sufhcc for (5.9) to hold. Applying (5.9) to gi{CnW) = E:^2_^(x^''-i) p 
C^'^) n W) and g2iCnW) = (A[,'''''^)2jy2(c('^-) n W), where C^^) denotes the union of k- 
faces of the d-polytope C, we see that (5.9) also holds for g — gi ~ g2- Thus, (5.8) is 
proved. From (5.8) and Lebesgue's dominated convergence theorem, we conclude that 



\fi'\s,t)\<E 



exp 



To show that fe^\s,t) becomes arbitrarily small as Wg grows large, we start with the 
obvious estimate 



2\W, 



exp 



0. (5.10) 



\fi'Hs,t)\<: 



( 2 

Ex exp{ist/f ) } - exp I - ^ Varx (4'^ {W,)) 



(5.11) 
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for fc = 1, . . . , d — 1. Next, we express C/g'^^ in terms of the centered measures ol^^\-) = 
Vk-i{Xn^~^^ n sl*^^ n (•)) — X'^''''^ VkC^Ti^ n (•)), which are conditionally independent given 
the initial tessellation X = {En}n>i- We have 

' ^' n>l 

and introduce, for fixed S >0, the conditional Lindeberg function 

' ^' n>l 

Further, for any £ > and 5 > 0, we define the events 

G,(e,5) = {iW(<5)<£} and H,ie) = {\LfHo) - (4'^^''^^ < e}. 

Since Lf\o) = YaTxizlf\Wg))/\We\, it follows from (5.8) that P(i?^(e)) — >g-,ooO- 
Below, we also need that P(Gg(e,J)) — >g_,ooO, following from the stronger result 

Lg{S) ^^g^oc 0, which we can show in the following way. Replacing (Wg) in (5.12) by 

n Wg) = i^k-iixi'^^^^ n Ei!"^ n Wg) + X^^''^^^ i^ki^it'^ n Wg) leads to the inequality 

Li'\a\Wgr/')<-^Y.Mw,nE„m^x{^lt\El^^ 

' ^ ' n>l 

for any a > 0. The set function g{C n W^) = Ex(Co''^(C('=) n W^))^l{5('=)(c(fe)nw)>a} 
translation invariant (due to the stationarity of Xq) and increases whenever W expands. 
Hence, g{-) fulfills the conditions needed to establish the almost sure convergence in (5.9). 
This implies that 

¥(^lhjijnpLgia\Wg\-'^')<^J^Ex{^l!'\c^'^^^^^ = 1. 

Consequently, by (3.5), P(limsupg^o^ Lg'^\S) < e) = 1 for any e > 0, that is, L^g\5) ^^g~>oo 
A suitable upper bound of the right-hand side of (5.11) can be obtained when both events 
Gg{e,5) and Hg{e) occur. From (5.11), it is easily seen that \fg^\s,t)\ does not exceed 
the sum 



El{G,(e,5)nff,(e)} 



Exexp{is;7f)}-cxp|-yL,(0) 



+ 2nG'gie,S)umis)). (5.13) 



We proceed with the factorization of the conditional characteristic function of Ug'^'^ 
given X, using the conditional independence of the random variables 9l!'\Wg), and 
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obtain 



Exexp{isC/f-)} = W ¥.xcM^s\We\-^'''e[^\w,)}. 

n>l 

Expressing the first equality in (5.8) by the centered measures 0^n\-), we have 
expj-^ Varx(4'Hw^.))} = H '^^pj^^'^^^^^-''^^^))' 

By means of the elementary inequality \xi ■ ■ ■ Xn — Vi ■ ■ ■ yn \ < — yi\ + ■ ■■ 
for complex numbers Xi,yi lying on the unit disc, we arrive at the estimate 



Ejcexp{isC/^''^^} -exp 



< 



{H„nVVe#0} 



'1" vm / 



Ex exp'^ IS — ^^^^^ (■ — exp 



(5.14) 



2\W,\ 



Further, using the well-known inequality |e^^' — X]fe=o ^^tl I — ^^nT ('^it^^ x e M^) for n = 2 



and n = 3, we find that, for any S > 0, 

( r. e^::\w,) \ 

- exp<^ IS — , > - 1 
V I J 



© + lT^Ex(e^(M^,))^ 



ell '^{\9\;:''(w,)\>sJ\w^\} 



(5.15) 



+ 



Analogously, applying the inequality |e^^ — 1 + a;| < for x > gives 



exp 



< 



(5.16) 



where we have used, in addition, that, for any n > 1 and ^ > 0, 
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Finally, combining the above estimates (5.14), (5.15) and (5.16) and taking into account 
both the abbreviation (5.12) and the fact that E,xdi!'\Wg) = 0, we find that 

Ex cxp{isC/f-) } - cxpj - (0) 

Regarding the latter inequality on the event {Gg{e,5) n Hg{e)}, we obtain from (5.13) 
that 

lim sup l/f ) (., t) I < .^e + + ^ {5' + e)) {e + {4'''^)') 

for arbitrary e,S >0. Thus, limg^oo fg^\s,t) = which completes the proof of (5.6). 
The proof of Theorem 5.1 ends with an application of (5.6) and the continuous mapping 
theorem (cf. [17]) to the linear combination (5.5), which proves (5.7). □ 



< .s^L W iS) + \fL^Lf^ (0) + ^L(^-) (0) (-5^ +Lf^ (6) ) . 



6. Examples of weakly dependent random 
tessellations 

There are only a few papers (e.g., [3, 10]) concerning weak dependence properties of 
stationary random tessellations apart from ergodicity. In fact, the assumption of ergod- 
icity turns out to be the weakest form of asymptotic independence of distant parts of a 
stationary tessellation X = {E!}„>i. Due to the individual spatial ergodic theorem (cf. 
[6]) ergodicity guarantees strong consistency for a number of intensity estimators based 
on a single observation in an expanding sampling window. To establish asymptotic nor- 
mality of these estimators, the distribution of X must satisfy certain mixing conditions 
expressed in terms of corresponding mixing coefficients. 

In the context of random tessellations X = {S}„>i in R'*, the a- and /3-mixing condi- 
tion have proved meaningful with mixing coefficients defined by 

aiAxiFi),AxiF2))= sup |P(yli n A2) - P(v4i)P(^2)|, 

AieAxiFi),A2eAxiF2) 

(6.1) 

PiAx{Fi),AxiF2))=E sup \P{A2\AxiFi)) -P{A2)\, 

A2eAx(F2) 

where Fi,F2 are disjoint closed subsets of R'' and Ax{F) denotes the cr-algebra gen- 
erated by the random closed set (Un>i^'^n n in the sense of Matheron (cf. [20] and 
also [10]). It is easily verified that a{Ax{Fi),AxiF2)) < l3iAxiFi),AxiF2)). However, 
the behaviour of both mixing coefficients is nearly the same for most of the models in 
stochastic geometry when the distance between Fi and F2 becomes large (cf. [12]). To 
verify Condition G, we are faced with two problems. First, to find, from the model as- 
sumptions, sharp bounds on the above mixing coefficients for Fi =F{a) := [—a, a]"* and 
F2 ■— G{b) —M.'^\ [—b,bY for b> a and, second, to prove a suitable CLT (or to use a 
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known CLT) for weakly dependent random fields whose assumptions follow from the 
derived estimates. For more details on CLTs of random fields and mixing conditions, 
including the influence of the dimension, the reader is referred to [23]. 

In [10], the /3-mixing coefficient l3{Ax{F{a))^Ax{G{a + r))) could be estimated for 
Voronoi tessellations in terms of the /3-mixing coefficient and certain void probabili- 
ties of the generating stationary point process of nuclei. In the special case of Voronoi 
tessellations generated by Poisson cluster processes with cluster radius i?o satisfying 
Eexp{/ii?o} < oo for some ft, > 0, the general bound decays exponentially in r. More 
precisely, it can be shown that 



r \ la 



f3iAxiF{a)),AxiGia + r)))<cA[-] + - cxp{-C2r} (6.2) 



for any r > 1 and a > 1/2, where the positive constants ci, C2 depend only on the dimen- 
sion d, h> and the intensity of the Poisson process of cluster centers. An estimate simi- 
lar to (6.2) holds for Poisson soft-core processes (cf. [29]), provided the soft-core radius Rq 
possesses an exponential moment. Furthermore, Condition G could be verified in [10] for 
stationary random tessellations X = {S}„>i in R'' satisfying E{h'k{X^''^ f] [0, 1)''))^+'' < oo 
for 1 < k < d — 1 and some S > and 

/3(^x(F(a)),^x(G(a + r)))<a'^-i/3i(r)l[i.,,](r)-f /32(r)l(,a,oo)(r) 

for any a > 1/2 and r > 1, where c > 2 is a constant independent of both a and r. 
Furthermore, /3i(-) and /32(-) are non-increasing functions on [l,oo) such that 



^2d. 



-i/3i(r) and V r^^'^ (/32(r))*/('+*) < oo. 

n — inn ' ^ 



Hence, from (6.2), we obtain, for a Voronoi tessellation X — {.=,}„>! in R'' generated by 
a stationary Poisson process with intensity 7 > 0, that 

\W,\-^l\vk{X^^'>^W,)~Tif{v)^^''-^^l''\W,\) ^ AA(0,rL(^)7^'-''^/') (6.3) 

n— i-oo 

for fc=l,...,(i— 1, where the mean value -Jif'iv) = E(z^fe(X('=) n [0, 1)'') and the asymp- 
totic variance f\JyV) refer to a PVT with intensity 7 = 1. The scaling rates in (6.3) 
are easily seen from the scaling property of the stationary Poisson process in R'' giving 
VkiX^''^ n Wg) = -f-^/'^Vkix''''^ n Wg^i/d), where X^''^ denotes the union of fc-facets of a 

PVT with unit intensity. An explict formula for the intensity ]I^^\i>) was first found (by 
R. Miles) to be 

^ (27r)'^-fe+lr(d-fc + fc/d) Hd(^d-k)+k-2 Kk-l ''"'^ 

{d-k + iy.d (Kd)*=/^ 

and derived in a different way in [21], page 64, whereas, for Tf,^{i'), no analytic ex- 

(2) 

pression is currently known. In the planar case, it is well known that ]l\ (i^) = 2 and the 
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^(^)-4A^^)^/A and = (r^^'^))^ + 1.6934(A^'V- 



approximate value 2(1^) = 1.0445685 was found in [1] by numerical evaluation of rather 
involved multiple integrals. 

We are now in a position to establish the CLT in (5.7) in a more explict form for the 
case of a planar X/Xo-ncsting with initial tessellation X being a PVT with cell intensity 
7 > and component tessellation Xq being either a PLT generated by a motion-invariant 
Poisson line process with intensity A > or another PVT with cell intensity A > 0. In 
both cases, we have 

„Z^-^(0' ^1.2)' (6-4) 

where 

From (3.2) we get a[,^'^' = cf 'aJ,^'^^ with cf = 2/7r and aJ,^'^^ = 2\/A if Xq is a PVT 
(cf. [24]), page 314 and a[,^'^' = A if is a PLT (cf (2.8)). To calculate 

{T^'^'^f=ll Var(.o(X«nC(i'))P"(dC), 

we first consider the case when Xq is a PLT with intensity A. Then, vq{xI^^ r^C^^'>) equals 
twice the number N{C) of Poisson lines hitting the polygon C. It is well known that 
N{C) is Poisson distributed with mean (and variance) XP{C)/tt, where P{C) denotes 
the perimeter of C. Hence, by EP(S*) === 47-^/2 (cf. [24], page 314) we obtain 

TT TT 

If Xq is a PVT with cell intensity A, we may again exploit the scaling properties of PVTs, 
giving 



jXEYarj^iiyaiXy HE " ')), 



where X and Xq are independent planar PVTs, both with unit cell intensity, and where 

11 of " 



(1) —*,(!; 

Xq n S denotes the finite set of points on the boundary of the typical cell of X 



induced by the 1-facets of Xq. A large-scale simulation study yields E Var^(i'o(-''^o ^ 
' )) = 2.7023. Summarizing the above results, we obtain the following expressions for 
yy^^'' and 2 (6-4), namely. 

' = -^7-^ and ^ ^ — V7^ + 1.6934A^ (if Xq is a PLT with intensity A), 

TT TT 

^ = — and rf 2 = 2.7023a/7A 6.7736A (if Xq is a PVT with intensity A). 

TT ' 



CLTs for functionals on facets of random tessellations 
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To conclude, it should be mentioned that Condition G can also be verified for a large 
class of Laguerre tessellations generated by Poisson-based point processes. The values 
of the variances in the previous formulae for higher dimensions can only be obtained by 
extensive simulation studies. Several generalizations of Theorems 4.1 and 5.1 are possible. 
For example, the manifold process X^f'^'' of {d— l)-facets can be replaced by the union 
of fc-facets of Xq for 1 < k < d — 2. In case of anisotropic component tessellations, the 
rose of directions of Xq is needed to express the mean and variance of zlf\Wg). 
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